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1 .  PROBLEM 


Optical  pulses  that  result  from  the  reflection  of  GaAs  laser  pulses 
from  water  clouds  have  been  measured.  The  amplitude  of  tlte  return  pulse 
as  a  function  of  time,  V(t),  depends  on  the  transmitted  laser  pulse, 
p(t);  the  overall  system  sensitivity  as  a  function  of  distance  (range 
response)  R(x)>  and  the  optical  properties  of  the  cloud.  For  the 

present  purpose,  the  properties  of  the  cloud  along  the  pencil  beam  path 

of  the  optical  probe  can  be  described  by  a  single  function,  C(x),  the 
cloud  signature. * Light  arriving  at  the  detector  at  time  t  must  have 
been  transmitted  at  an  earlier  time,  t  -  t  (t  >  0).  t  is  the  time  that 
it  takes  light  at  speed  c  to  travel  from  the  transmitter  to  the  point  of 
reflection  (distance  *  x)  and  then  back  to  the  detector,  which  is  near 
the  transmitter.  Therefore, 

t  =*  2x/  c-  >  ( 1 ) 

All  pairs  of  t  and  x  such  that  t  0  and  equation  ( 1 )  holds  will 

contribute  to  the  return  at  time  c ,  with  the  contribution  weighted 
according  to  the  system  range  response  at  the  distance  x  and  the  cloud 
signature  at  that  location.  Therefore,  the  total  return  is  given  by  the 
convolution 

V(t)  =  K  /g  P(t  -  t)R(x)C(x)  dx  ,  x  «  cr/2  ,  (2) 

where  K  is  for  normalization. 

The  purpose  of  this  study  is  to  investigate  a  method  for  finding 
C(x)  given  V(t),  P(t),  and  R(x).  P(t)  and  R(x)  are  constant  functions 

of  the  system.  V(t)  has  been  measured  for  a  large  number  of  pulses. 

One  method  for  solving  for  C(x)  in  equation  (2)  is  as  follows.2 
Using  the  notation  of  McGuire,2  let 

h(t)  *  K* C( ct/2 ) • R( ct/2 )  .  (3) 

Now,  to  consider  equation  (2)  numerically,  let 

t  =  n*At  ,  r  =  i*At  • 


H .  Burroughs ,  Computation  of  Cloud  Backscatter  Power  as  a  Function 
of  Time  for  an  Active  Optical  Radar  (U ) ,  Naval  Weapons  Center  NWC  TP 
5 090  (April  1971).  (CONFIDENTIAL) 

2 Dennis  McGuire  and  Michael  Conner,  The  Deconvolution  of  Aerosol 
Backscattered  Optical  Pulses  to  Obtain  System- Independent  Aerosol 
Signatures ,  Harry  Diamond  Laboratories  HDL-TR-1944  (June  1981). 


The  integral  then  becomes  a  buru 

n 

V(n*At)  -  At  *<<»»  “  i)At)  <h(  iAt)  *  (4) 

The  upper  limit  equals  n  because  P  *  0  for  t  <  0.  Thinking  of  V  and  h 
as  vectors ,  the  notation  becomes  “* 

$  -  At»pft  ,  (5) 


The  problem  of  solving  for  C(x)  (by  here  solving  for  h)  then  becomes  one 
of  matrix  inversion.  The  solution  to  equation  (5)  is 

ft  -  (1/At)p“lV  .  (6) 

Equation  (6)  does  not  give  satisfactory  results^  if  V  contains  noise. 
This  failure  is  probably  because  the  form  of  P(t)  causes  some  large 
elements  to  appear  in  P"1 ,  giving  too  much  significance  to  same  small 
variation  in  V  due  to  noise.  It  is  not  surprising  that  there  is  some 
difficulty  related  to  dealing  with  P”1.  This  difficulty  is  because  the 
determinant  of  P  is  [P(At)3n»  which  is  very  small  since  P(At)  is  the 
first  nonzero  point  of  the  transmitted  pulse.  One  might  easily  expect 
this  problem  from  another  point  of  view;  namely,  since  the  convolution 
of  C  (eq.  2)  will  smooth  out  small  bumps  in  C,  the  deconvolution  of  V 
using  the  same  equation  will  badly  exaggerate  small  bumps  (noise). 


2Dennis  McGuire  and  Michael  Conner,  The  Deconvolution  of  Aerosol 
Backscattered  Optical  Pulses  to  Obtain  System-Independent  Aerosol 
Signatures,  Harry  Diamond  Laboratories  HDL-TR-1944  (June  1981), 


A  way  to  avoid  all  such  difficulty  ia  to  Mk*  tha  approximation 
that  P(t)  *  A*S(t  -  tQ) ,  whara  A  ia  a  normalisation  constant  and  t^  is 
tha  cantar  of  tha  transmitted  pulse.  Than  station  (2)  immediately 
yialda  r/  v-> 


This  discards  all  information  about  tha  ahapa  of  tha  tranamittad  pul an, 
thus  reducing  tha  rasolution  in  C(  x)  to  about  1  to  2  m  ainca  tha  trana¬ 
mittad  pulaa  has  a  full  width  at  half  maximum  ( FWHM)  of  about  6  ns. 


2.  DESCRIPTION  OF  METHOD 

In  this  investigation ,  aquation  (?)  is  usad  as  a  first  approximation 
to  C(x).  Since  one  desires  a  cloud  signature  with  batter  than  tha  1-  to 
2-m  rasolution  achieved  by  using  equation  (7),  soma  effort  is  than  made 
to  modify  C(x)  to  find  a  more  accurate  cloud  signature*  In  this  discus¬ 
sion,  a  more  accurate  (or  "better")  cloud  signature  refers  to  a  signa¬ 
ture  with  a  smaller  error,  where  error  is  defined  as  follows*  Insert 
the  cloud  signature  currently  being  considered  into  equation  (2)  to 
calculate  the  return  pulse  (Vc,  subscript  c  for  calculated)  that  would 
result  from  that  signature.  Oompare  that  with  the  measured  return  pulse 
(Vm)  defining 

error  « 

Various  functions  (as  discussed  later)  are  tried  for  C(x),  and  the  C'x)' 
that  gives  the  smallest  error  is  recorded  as  the  extracted  cloud  signa¬ 
ture. 


Now,  equation  (6)  immediately  gives  the  cloud  signature  with  error  = 
0,  but  that  signature  is  noise  dominated  nonsense.2  Here,  then,  one  is 
not  seeking  the  absolute  minimum  in  the  error  (which  would  be  zero) ,  but 
rather  one'  seeks  a  relative  mininum  in  error  by  varying  C(x)  in  seme 
gentle  way  about  the  6  (delta)  extracted  (eq.  7)  first  guess. 

The  following  measures  are  available  to  prevent  the  cloud  signature 
from  becoming  wildly  bumpy: 

a*  Restrict  how  far  the  successive  estimates  of  C(x)  can  vary  from 
the  original  6  extracted  C(x). 


2 Dennis  McGuire  and  Michael  Conner,  The  Deconvolution  of  Aerosol 
Backscattered  Optical  Pulses  to  Obtain  System-Independent  Aerosol 
Signatures ,  Harry  Diarth-jnd  Laboratories  HDL-TR-1944  ( June  1981), 


b.  Subject  the  return  pulse ,  Vm(t),  to  e  low-pass  filter. 

c.  Subject  the  final  answer  for  C(x)  to  a  low-pass  filter. 

The  severity  of  each  of  these  restrictions  Is  easily  varied.  It  was 
hoped  that  with  theae  restrictions  a  relative  minimvm  In  the  error  could 
be  found  associated  with  some  reasonable  cloud  signature  that  would  have 
the  same  gross  properties  as  the  delta  extracted  signature,  but  with 
sharper  features.  Before  sample  calculations  were  carried  out.  it  was 
not  clear  whither  restrictions  and  filtering  (measures  a  to  c)  adequate 
to  smooth  out  the  noise  would  simultaneously  doom  resolution  to  worse 
than  1**  to  2^  resolution  of  the  delta  '  function  case.  If  filtering 
alone  (b  and  c)  can  successfully  suppress  the  noise  in  a  given  data  bank 
while  still  allowing  improved  resolution  (this  would  depend  chiefly  on 
the  signal  to  noise  ratio  for  that  data),  then  one  may  consider  using 
equation  (6)  after  all,  with  the  appropriate  filtering  of  V(t)  and 
C(x).  On  the  other  hand,  if  the  a  priori  assumption  that  the  desired 
cloud  signature  closely  resembles  the  delta  extracted  signature  is 
essential,  then  the  msthod  described  below  msy  be  useful. 

In  the  present  method,  the  calculations  are  done  by  computer.  The 
cloud  signature,  C(x),  is  in  digitised  form  with  100  values  for  C. 
These  values  are  at  0.1 5-m  intervals  in  distance  (x),  covering  a  total 
of  15  m.  The  first  estimate  of  C(x)  is  provided  by  the  delta  function 
method;  that  i9,  equation  (7)  is  used  for  100  valxies  of  t  (corresponding 
to  100  values  of  x  with  x  -  ct/2).  The  various  values  of  C(x)  are  then 
varied  in  an  effort  to  find  a  C(x)  with  a  smaller  error  (as  defined 
above).  The  variation  of  C  proceeds  as  follows;  For  a  chosen  initial 
value  of  Q  (Q  >  1),  begin  with  the  first  point  )  -  0.15  m]  and 

consider  the  following  five  possibilities: 

a.  Leave  C(x)  unaltered. 

b.  Multiply  by  Q. 

c.  Divide  C(x^)  by  Q. 

d.  Multiply  C(x.)  by  Q  and  decrease  the  next  value,  C(x^+^) ,  by  the 
amount  by  which  C(x^)  increased. 

e.  Divide  C(x±)  by  Q  and  increase  C(xi+1)  by  the  amount  by  which 
Cfx^)  decreased. 

(Choices  d  and  e  are  motivated  by  the  consideration  that  they  leave  the 
integrated  return  less  changed;  since  the  original  guess  had  approx* 
imately  the  correct  total  energy,  these  possibilities  are  probably 
desirable.)  One  then  takes  whichever  of  the  five  resulting  cloud  signa¬ 
tures  that  has  the  lowest  associated  error  and  considers  it  to  be  the 
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current  (or  "new  and  improved")  cloud  signature »  this  variational 
method  is  then  repeated  for  the  next  point  and  successive  points  through 
C(x100)  ■  15  m.  this  constitutes  one  pass.  If  after  any  pass  the  error 
is  smaller  than  it  was  at  the  beginning  of  th^  pass,  it  is  deemed  worth¬ 
while  to  make  another  pass  by  using  a  smaller  value  of  Q, 

°new  *  <n<3  +  O/tn  +  *>  ,  n  >  0  , 

to  achieve  finer  variations  in  C(x).  The  flow  chart  and  coding  for  the 
computer  program  that  carries  out  this  procedure  are  shown  in  appendices 
A  and  B,  respectively.  The  values  for  n  and  the  initial  Q  can  be 
specified  as  desired  for  each  computer  run.  The  particular  values  of 
( «  ^initial*  n  that  were  used  to  generate  the  examples  in  this 

report  are  given  in  appendix  C.  The  cloud  signature  is  constrained  to 
vary  only  within  a  certain  region  because  no  value  of  C  could  be 
multiplied  or  divided  by  more  than  approximately  Qi*Q2*Q3*...*S^,  and 
this  product  has  a  finite  value  depending  on  Qj  (»  Qinitial^  and  n  ^ae* 
app  D) . 

In  this  wsy,  C(x)  varies  until  some  relative  minimum  in  the  error  is 
found.  When  a  pass  is  executed  (one  Q,  all  100  points)  with  no  decrease 
in  the  error,  the  process  is  stopped  and  the  cloud  signature  is 
recorded. 

3.  RESULTS 


The  characteristics  of  the  signature  found  by  this  method  depend 
critically  on  the  signal  to  noise  ratio  of  the  return  signal,  V(t),  as 
the  following  examples  show. 

3. 1  Case  I 


An  idealized  cloud  with  the  backscatter  coefficient  propor¬ 
tional  to  the  extinction  coefficient3  and  extinction  coefficient  profile 
as  shown  in  figure  1  results  in  the  cloud  signature  shown  in  figure  2. 
Here  one  assumes 1  * 2 


C(x) 


-2/*0(s)  ds 
U(x)e  ° 


.  H,  Burroughs,  Computation  of  Cloud  Backscatter  Power  as  a 

Function  of  Time  for  an  Active  Optical  Radar  (U ) ,  Naval  Weapons  Center 
NWC  TP  5090  (April  1971 J.  ( CONFIDENTIAL ) 

2Dennis  McGuire  and  Michael  Conner,  The  Deconvolution  of  Aerosol 
Backscattered  Optical  Pulses  to  Obtain  System-Independent  Aerosol 
Signatures,  Harry  Diamond  Laboratories  HDL-TR-1944  (June  1981), 

3ZJ.  Diermendjian,  Electromagnetic  Scattering  on  Spherical  Poly¬ 
dispersions,  American  Elsevier,  New  York  (1969),  The  ratio  of 
backscatter  coefficient  to  extinction  coefficient  depends  on  the 
wavelength  of  the  radiation  and  the  aerosol  particle  size 

distribution .  If  the  latter  is  a  locally  homogeneous  property  of  the 
cloud,  y/o  (backscatter/extinction  coefficients)  will  be  constant. 
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The  return  pulse  (calculated  by  a  straightforward  convolution,  eq  2) 
that  would  result  from  such  a  cloud  was  subjected  to  the  present  signa¬ 
ture  extraction  algorithm.  The  algorithm  includes  filtering  out  high 
frequency  bumps  in  the  return  pulse  and  signature.  The  resulting 
extracted  signature  (fig.  3)  is  relatively  close  to  the  actual  signature 
and  provides  better  resolution  of  the  sharp  peak  than  the  delta  function 
extracted  signature  (fig.  4). 

Noise  in  the  return  signal  becomes  exaggerated  on  deconvo¬ 
lution.  This  exaggeration  means  that  the  signature  extraction  becomes 
less  reliable  as  the  signal  to  noise  ratio  decreases.  To  demonstrate 
this  effect,  the  return  pulse  from  the  idealized  cloud  has  been 
corrupted  with  noise  and  the  signature  extraction  has  been  attempted 
again.  The  results  for  a  signal  to  noise  ratio  of  about  100  to  1  are 
shown  in  figures  5  and  6  and  for  a  signal  to  noise  ratio  of  about  40  to 
1  are  shown  in  figures  7  and  8.  For  a  signal  to  noise  ratio  of  100  to 
1,  the  initial  rise  and  sharp  peak  are  just  slightly  more  accurate  than 
the  delta  function  was,  but  a  large  amount  of  noise  has  shown  up  in  the 
exponentially  decaying  tail.  For  a  signal  to  noise  ratio  of  about  40 
(fig.  6),  the  situation  is  worse. 


X  (m) 

Figure  1.  Case  I;  extinction  coefficient. 
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Zm.  Jii-f  fc;  r*nt 


a  (m 


Figure  8.  Case  Is  delta  extracted  cloud  signature 
(return  pulse  signal  to  noise  ratio  ~  40). 

3.2  Case  II 


Figures  9  to  16  follow  the  same  format  as  figures  1  to  8,  but 
for  a  different  idealized  cloud,  case  II.  Here  the  improvement  is  less 
substantial  because  the  case  II  cloud  has  a  softer  edge,  so  the  delta 
extracted  signature  is  closer  to  the  correct  signature. 


Figure  9.  Case  II:  extinction  coefficient 


Figure  11.  Case  II:  extracted  cloud  signature 
(noiseless  return  pulse). 
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figure  16.  Case  II:  delta  extracted  cloud  signature 
(return  pulse  signal  to  noise  ratio  ~  40). 


4.  CONCLUSIONS 

If  the  data  have  a  signal  to  noise  ratio  of  about  100  to  1  or 
better,  it  is  possible  to  achieve  slightly  better  resolution  in  the 
cloud  signature  than  the  resolution  that  would  be  available  from  equa¬ 
tion  (7),  dividing  by  the  range  response  (delta  function  extraction). 
The  extent  of  the  improvement  in  the  resolution  depends  on  the  signal  to 
noise  ratio.  There  is  no  evidence  that  restriction  a  isect.  2) 
( limiting  how  far  estimates  of  the  cloud  signature  can  vary  from  the 
delta  function  estimate)  is  essential  to  the  reasonable  behavior  of  the 
cloud  signature.  Low-pass  filtering  of  the  return  pulse  or  resulting 
cloud  signature  probably  is  necessary.  This  necessity  indicates  that, 
with  a  signal  to  noise  ratio  greater  than  100  to  1  and  adequate 
filtering,  perhaps  equation  (6)  should  be  reconsidered  as  a  method  for 
solving  this  problem.  It  may  be  worthwhilte  to  consider  using  a  hybrid 
method  by  which  equation  (6)  could  be  used  to  treat  that  portion  of  the 
return  where  the  signal  is  high,  while  settling  for  the  delta  extracted 
signature  where  the  signal  is  lower. 


It  was  considered  desirable  to  analyse  a  select  portion  of  the 
measured  data  by  using  the  present  method.  Because  of  a  computer  mal¬ 
function,  thJ s  analysis  has  not  been  done.  It  is  suggested  that  for 
each  return  pulse  two  estimates  of  the  cloud  signature  be  recorded,  one 
using  the  algorithm  discussed  in  this  report  and,  for  comparison,  one 
using  the  delta  function  extraction.  Unlike  the  test  cases  shown  in  the 
figures,  th^  real  data  provide  no  answer  key  and,  hence,  it  may  not  be 
clear  which  signature  is  more  accurate.  However,  it  may  prove  instruc- 
ti ■.’«  to  generate  graphs  of  the  two  estimates  for  various  pulses  and  to 
not®  the  similarities  and  differences  in  the  two. 

Interestingly,  McGuire2  predicts  that  the  signal  to  noise  ratio 
should  decrease  by  a  factor  of  about  0.14  on  deconvolution,  with 
filtering  of  the  return  pulse  as  was  done  in  the  test  cases.  This 
compares  the  signal  to  noise  ratio  of  the  filtered  return  with  that  of 
the  cloud  signature.  Tn  the  present  study,  the  signal  to  noise  ratio  in 
the  unfiltered  return  signal  is  compared  with  that  of  the  cloud  sig¬ 
nature.  Also,  the  cloud  signature  was  smoothed  here.  For  these  two 
reasons,  one  expects  the  apparent  increase  in  noise  to  be  less  than  that 
predicted  by  McGuire.2  The  tvro  test  cases  here  show  an  increase  in 
noise  (compared  with  signal)  on  deconvolution  by  a  factor  of  about 
three,  roughly  consistent  with  this  expectation. 
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APPENDIX  A. --PLOW  CHART  OF  COMPUTER  PROGRAM  FOR  CLOUD 
SIGNATURE  EXTRACTION 

A  computer  program  was  developed  to  deconvolve  an  optical  return 
pulse  in  such  a  way  as  to  minimise  the  effect  of  return  signal  noise. 
The  flow  chart  for  this  computer  program  follows  (fig.  A-1). 


Figure  A-1.  Flow  chart  of  computer  program 


APPENDIX  8. — COMPUTER  CODS  FOR  CLOUD  SIGNATURE  EXTRACTION 

A  computer  program  was  developed  to  deconvolve  an  optical  return 
pulse  In  such  a  way  as  to  minimise  the  effect  of  return  signal  noise.  A 
FORTRAN  listing  of  this  program  follows. 
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C  CLOUD  SIGNATURE  IHTMCTIOH 

C 

C 

C 

C 

C 
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UC<M. Z  ) •  RETURN  PULSE  AT  EACH  NSEC  THAT  WOULD  RESULT  FROM  THIS  SIGNATURE. 
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BMUBBM*CMNOM 

4M 

00  TO  1BB0 

4H 

c 

437 

e 

BLOW  CHMMT  BOM  MB. 7 

4M 

B 

4M 

3B7B 

BO  30TB  X-l.M 
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fl 

M 


441 

44® 

443 


44® 

447 

449 

443 

4lt 

4«1 


e 

e 

e 

e 


4«3 

4«4 

411 

4U 

4*7 

411 

4*3 

4*3 

4*1 

4tl 

4*3 

4*4 

4C* 

4SC 

4*7 

4*3 

4*9 

47* 

47t 

47* 

473 

474 
47* 
47* 
477 
47* 
479 
4*9 
4*1 
4«« 
4*3 
4*4 
4*9 
4*« 

**n  > 
4*7 

4*9 

4*9 

49* 

491 

9CTU9H' 

49* 

493 


VU9d*4>-0d#l> 

CONTtNg* 

KIWI 

I-UNIN9**-9NIPT*3  TIN* 

*  90  3974  1-11*9* 

41*1 

T*NP 1  -9 1  -9 .  **T  into 
T*N49-P 1  .  «T1I*0 

Viiu-Ttnoi 

ip  tiim.uT .  t  )Uiu  •  i 

UWTHPI 

lPtUl"9.0T.**>U«*-9* 

Ctl.l >-*• 

^So'Sir*  K-um.una 

T*H9-4*#dPl-ro/*TiH«o> 

CU.I  l-C«l.l»4VUm»C.4>»tK9«-71l9*> 
PN04PI  -PN04N-CXP  <  -TINP  ) 

»"  JSCI'SS.t.tv—w 

3974  CONTXNU* 


0 

c 


PLOU  CH49T  90M  40.9 
*un-a. 

DO  3*99  1-1*19*  _  . 

*un-*u«4wn<x>*wmu 
CONTINUE 

*44094-*9*09/9un 

D£LC«A-t)ttt4N'*yn _ _ 

XP< *94094. LT.EC4XT >90  TO 

PLOU  CHANT  90N  4* 

C  1CLACS-3 

XCT93-XCT9341 
3199  CONTINUE  ^ 

IP«14LOTP.N*.1>00  TO  39*0 
00  3193  Z-1.9* 

J-1M9 

CSHlPTil  >-Cd*J> 

CALL  fU9f»LT(X9HXPT*99* ‘  <',1»C9H1PT» ,*19*,3«1P9 

319*  CONTINUE  ^ 

DO  3197  1-1*100.. 

OCPLOTd  )*UC<  1 . 1 ) 

,W  c2lL1£5pLT(TN4TW<.14*.  ^C*.4,WC4L0T.»C4tC. 

*19. XP44N1 ) 

*U*  ^WiTcTNATj:.!**.-^-. 

C4.UR. *HK49»  9*71901* •  19* 1P94H1  ) 


A 


ikjMUil*  ....  "i  ifmiii <f  liMilrttliXl 
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MM  CONTINUE 


caaaaaaaaaa 

cminiiit 

caaaaaaataa 

caaaaaaaaaa 


Jin  wixrriNi 

MERE.  UAXTt  TO  TAPE*  0(1,X>,  1*XX,M 
ITHXS  It  THE  MODIFIED  UCRSIOM  OP  T>« 
CLOUD  SXONATURE.) 


OC  TO 


end  or  procrah  Dunmutv 


CALL  SCREEN 

UAITC<  XOUTDW.MMl  JIHAHV 
URXTE<  XOUTDU.MSMlLOUCTR 
UAXTE<  XOUTDU. 90304  >  ICTR1 
UOXTC<  X  OUT  CO,  MMX )  XCTOt 
UA2TC<  XOUTDU.  MM?  )XCYM 


esastaattas 

caaaaaasaaa 

caaaaaaataa  XNMRT  A  CALL  rm  MERE. 

caaaaaaaaat 

caaaaaaataa 

e 

c 

c 

stop 

AM  XCTR1-XCTRU1 
CO  TO  300 

AMO  UAXTE<XOUTDU»MWO>XTERP 
STOP 

c 

C  xr  UE  OET  A  SAD  FRAME  I 

c 

SMO  UAXTEt  XOUTDU.MSSl  JXSTAT 
C 

c 

Ctaaaaaaaaa 

caaaaaaaaaa 

caaassaaaaa  INSERT  A  CALL  FFHN  HERE. 

caaaaaaaaaa 

caaaaaaaaaa 

c 

0 

e 

c 


c 

c 
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Mill  FOMPT!'  ENTER  CRITIOAL  WHOM.*) 

Mill  F0RMAT!E1B.3> 

Min  format!  •  enter  initial  value  for  o.  ’  > 

Mill  FORMAT!'  ENTER  N  (UCIQHT  FAC  TON  NON  0  DECREASE  >. '  1 
NIM  FORMAT!  '  CNTKN  OUTPUT  TANS  DEVICE.') 

Mill  FORMAT!'  ENTER  INTEGRATED  PULSE  UtlOMT  KLOU  UNION 
C'SIONAL  SHOULD  NOT  BE  PROCESSED. *  > 

Milt  FORMAT C  *  ENTER  NUMBER  OP  PULSED  TO  BE  PROCESSED. ' > 

■Mil  FORMAT (  '  PON  FILTERING  OP  RETURN  RULES.  ENTER  ' 

C.'l'E  TIME  (IN  NSEC  >*  '/‘DO  NOT  USE  ZERO*/ 

C' SHOULD  USE  AT  LEAST  .f  NSEC  SO  YOU  AUQ.'/ 

C'OUER  MORE  THAN  ONE  RASTER  POINT.'} 

Mill  FORMAT ( '  POR  FILTERING  OP  CLOUD  SIOMATURE.  ENTER  ' 

C.'l'E  DISTANCE  (IN  UNITS  OP  .IB  METERS)*  '."DO  NOT  UBC  KRO') 
SOSES  FORMAT «'  ENTER  X.  FOR  PLOTS.') 

80301  FORMAT! IX.  IS.  '  RETURN  PULSES  UCRE  EXAMINED. '  ) 

00303  FORMATdX.lt.'  OF  THESE  PULSES  UERE  BELOU  THRESHOLD.') 

80304  FORMAT! IX. IS. '  OF  THESE  PULSES  HAS  LESS  THAN  ISO  NSEC'/ 

C'OF  USEABLE  RETURN.') 

SOSOS  FORMAT! IX. IS. ‘  CLOUD  SICNATURES  UERE  SUCCESSFULLY  EXTRACTED.') 
B038T  FORMAT! IX. IS. '  CASES  FAILED.') 

SSSSS  FORMAT!'  USE  DECIMAL  POINT') 

8SSSS  FORMAT (F IS. 3 ) 

SOSOS  FORMAT C9AS) 

90US0  FORMAT!'  ASSIGNMENT  ERROR '.10. 'FROM  CETFR' ) 

00SS1  FORMAT!'  IF  CODE  »  130.  FILE  MARKER  FOUND.'/ 

C'  ANY  OTHER  CODE  MEANS  TAPE  ERROR.'/ 

C*  CODE  •  '.110) 

END 

SUBROUTINE  SUBPLT ! PASS1 . IPASSS. IPASS3. IPASS4. 

CPASSS. IPASSS. I PASS?. IFPAME ) 

DIMENSION  PASS1 ( X ). IPASSS! 4 ). IPASSS! 4 ). PASS*! X ) 

C 

C  VARIABLES  PASSED* 


1  «  X  ARRAY  NAME 

a*  •  points 

3*  'X  LABEL' 

4*  •  CHARACTERS  in  label 

Et  V  ARRAY  NAME 
S*  'V  LABEL' 

?•  t  CHARACTERS  IN  LABEL 

IFRAMEt  FRAME  CODE 

CALL  SCREEN 

URITE! 1 . 90S ) XFRAME 

CALL  SCALE ( PASS! . IPASSS. 0 . .1.0) 

CALL  SCALE(RASSS.IPAS08.0.*1.1) 

CALL  ENTQRA 

CALL  XAXXf! IPASSS. IPA804.8. ) 

CALL  VAXIS! IPASSS. IPASST.S. > 

CALL  DATAQ ( PASO). PASSS, IPASSS, 1.1) 
CALL  CXI TOR 
READ (1.903) I FAKE 
RETURN 

•03  FORMAT!  ID 

MB  FORM  AT  (SIX.'  FRAME'.  IB) 

END 
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The  return  pulse  was  filtered  by  replacing  each  digitised  value, 
V(t^),  by  a  weighted  average  of  the  points  in  its  neighborhood.  The 
weight  factor  decreases  exponentially  as  t  moves  away  from  t,.  The  time 
(in  nanoseconds)  by  which  the  weight  factor  has  decreased  to  1/e  is 
called  t0«  The  filtering  of  the  cloud  signature  is  similar,  except  that 
the  units  for  xQ  are  0.15  m  (corresponding  to  1  ns). 

The  filtering  values  and  the  values  for  and  n  used  to  generate 
the  figures  in  the  main  body  of  this  report  were  these) 

Case  I 

tQ  “  1.3  (”  1.3  ns)  , 
x0  «  1.2  («  0.195  m) , 

Qi  *  U2t 
n  ■  2. 

Case  II 

tQ  1.3  (»  1.3  ns)  , 
xQ  =  1.3  («  0.195  m) , 

-  1.1, 
n  =2. 
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As  discussed  in  the  main  body  of  this  report,  the  extent  of  vari¬ 
ation  of  the  cloud  signature  depends  on  the  product  of  the  Q^'s.  The 
author  has  not  been  able  to  write  this  product  in  terms  of  and  n  in 
closed  form,  but  some  discussion  and  numerical  results  are  given  below. 

Notes  1  <  Q  <  2,  n  >  0,  0  <  Z  <  1,  0<a<  1. 


Find 


n  Qi  =  Ql*Q2*Q3*Q4  •  •  .*  » 


i=1 


where 


Let 


1  +  nQi 
3i+1  “  n  +  1 


zi  =  Qi  "  1 

Then  from  equation  (D-2)  it  follows  that 


i+1 


■  rf-r  (*i)  -  • 


n  +  1 

Using  equation  (D-4)  repeatedly  shows  that 

zi  =  ai”1zi 

Back  to'  the  task  at  hand, 


n  %  -  n  k  ♦ »  . 

i=1  1  i=1 

ln (noi)  =  in  [n(zi  + 1)]  - 1>  (z±  + 1) 


( D-1 ) 


(D-2) 

(D-3) 

(D-4) 
(D-5 ) 

(D— 6 ) 


(D— 7 ) 

(D-8 ) 


In  (riQi)  =  In  (l  +  Zx]  +  In  (l  +  oZx)  +  ln  (l  +  a2Zi) 


(D-9) 


+  In  (1  +  a3Zi)  +  .  .  ,  . 


Using 

In  (1  +  x)  =  x  -  +  2£_  -  ^  ^  .  ,  .  for  x2  <  1  and  x  =  1 

in  each  term  in  equation  (D-9)  yields 
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In  C  n  Qt)  = 


r  Z?  Z?  ^ 

LZl  2  +  “3  -  “4  +  T  *  *  *_ 

r  m*  (azi)3  (°zi )-  w*  i 

+  jazi 2 +  —5 - 3 +  —5 •  • 

r,  (a)2  (a)3  (a)*  (a)s 

+  La2zi - r—  +  —3 - 4—  +  5 - 

[  (a3Zl)2  (a3Zl)3  (a3^  (a3Z  ^ 


u  zi - —  +  — T 


+  •  «  • 


Collecting  like  powers  of  Z, 


(nQj  “  Zi  +  a  +  a2  +  a3  +  •  • 

-  \  [1  *  a2  «■  a-  +  06  *  •  •  -J 

+  T  £  +  0,3  +  “6  *  “5  +  •  •  •] 

-  [i  ♦  •"  «8  *  °12  ♦  •  •  •] 


+  •  •  • 


Summing  the  algebraic  series  in  each  square  bracket  yields 

t~j  p  -1  r  _ 

r^=] ' T  [r^J '  [rh^ 


^  •  • 
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These  terms  get  very  small  very  fast.  The  sum  therefore  can  be 
approximated  by  the  first  few  terms.  The  exponential  of  the  sum  then 
gives 

Jlt«i  • 

Table  D-1  was  generated  by  computer.  Using  the  first  five  terms  of 
equation  (D-12)  gives  the  same  results  as  the  first  100  terms  of  equa¬ 
tion  (D-1)  to  the  accuracy  shown  here. 


TABLE  D-1.  NUMERICAL  RESULTS  TOR  PRODUCT  FI  Q, 

i-1 


1.4 

1.3 

1.2 

1.1 

1.05 

1.02 

1.01 

i  i 

1  3 

1.70 

1.50 

1.32 

1.16 

1.08 

1.03 

1.02 

-?  4 

i  1 

|  \ 

2.04 

1.73 

1.46 

1.21 

1.10 

1.04 

1.02 

2.95 

2.29 

1.76 

1.34 

1.16 

T.  06 

1.03 

-  4 

i  ] 

4.25 

3.04 

2.  14 

1.48 

1.22 

1.08 

1.04 

6.13 

4.02 

2.58 

1.63 

1.28 

1.10 

1.05 

'  1 

8.84 

5.32 

3.  13 

1.79 

1.34 

1.  13 

1.06 

t  1 

-i  Jj 

55.1 

21.6 

8. 12 

2.92 

1.72 

_ 

1.24 

1.12 
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